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Abstract. The dynamic linear response theory of a general Ising model weakly 
coupled to a heat bath is derived employing the quantum statistical theory of Mori, 
treating the Hamiltonian of the spin bath coupling as a perturbation, and applying 
the Markovian approximation. Both the dynamic susceptibility and the relaxation 
function are expressed in terms of the static susceptibility and the static internal 
field distribution function. For the special case of the SK spin glass this internal 
field distribution can be related to the solutions of the TAP equations in the entire 
temperature region. Application of this new relation and the use of numerical 
solutions of the modified TAP equations leads for finite but large systems to explicit 
results for the distribution function and for dynamic linear response functions. A 
detailed discussion is presented which includes finite-size effects. Due to the derived 
temperature dependence of the Onsager-Casimir coefficients a frequency-dependent 
shift of the cusp temperature of the real part of the dynamic susceptibility is found. 
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1. Introduction 

The Ising model of Sherrington-Kirkpatrick (SK) with quenched random bonds is 
the most important representative of a class of long-ranged models all describing spin 
glasses (for general references see [Q, || For the static analysis of this model two 
complementary but conceptually different approaches exist. The first approach uses 
the replica method |l| and the breaking of the replica symmetry || to calculate bond- 
averaged quantities. The approach of Thouless-Anderson-Palmer (TAP) |J is based on 
more conventional techniques and does not perform this bond average as it is expected 
that macroscopic physical quantities will be independent of the particular configuration 
in the thermodynamic limit. 

Although the TAP equations are well established 0, [| [7|] they are still a field of 
current interest. It is suspected that not all aspects of this approach have yet been 
worked out. Recently the author || has reanalyzed the stability of these equations 
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and has shown that unstable states cannot be described by the original TAP equations. 
Therefore he proposed modified TAP equations which turned out to be useful for explicit 
numerical calculations of the characteristic features of the SK model of finite-sizes ||. 

Dynamical questions are of great importance for the physics of spin glasses (for 
reviews see [0. |3|. ||). Therefore numerous dynamical extensions have been added to the 
SK model. Following the early dynamical approaches |T], [TI| Glauber dynamics has been 
used by various authors ||, |TTJ, pl| . In addition, Langevin dynamics has been employed 
in the studies [|l3j mostly for the soft spin version of the SK model. Both the Langevin 
and the Glauber dynamics are basically phenomenological and can at best be justified 
partially by microscopic arguments. Note this also applies to the work of Szamel [|12] 
although this approach is formulated in the 'spirit' of the microscopic Mori theory. 

It is the aim of this paper to present a compete microscopic analysis of the dynamical 
linear response for the SK model coupled to a bath. An adequate tool for this purpose 
is the general theory of Mori |L4|, [15| . This quantum-statistical approach is used in 
the present work and is worked out not only for the SK model but for a general Ising 
model. Such a treatment is obvious and straightforward. Nevertheless, to the authors 
knowledge it has previously not been published . 

As usual the results of the Mori theory are expressed in terms of static equilibrium 
quantities. Together with the static isothermal susceptibility, it is the internal field 
distribution function [|18| , |19f which fully determines the dynamical linear response. At 
this point we restrict the approach to the SK model. It is shown that the internal field 
distribution function can be related to the solutions of the TAP equations. Employing 
the approximate TAP solutions ||, all the quantities of the linear response theory can 
be numerically calculated in the entire temperature regime for all external fields. 



The references [ 15 , |T9| | showed that an exact and complete description of the 
thermodynamics of Ising models can be formulated in terms of the internal field 
distribution function. Thus as a byproduct the presented results of this function are of 
some interest independent of the dynamical linear response problems. 

Following a description of the microscopic Hamiltonian, the Mori approach for a 
general Ising system is performed in section 2. The internal magnetic field distribution 
function for the SK model is treated in section 3. Both the analytical and numerical 
results for the field distribution function, the dynamic susceptibility and the response 
function are explicitly presented in section 4. Finally, some concluding remarks can be 
found in section 5. 

2. Linear response for a general Ising system 

2.1. The microscopic system 



i 

2 

spins interact via an arbitrary Ising spin-spin interaction Jy(= Jji) and are described 



A system of iV spins Sj with s = \ is considered in the presence of external fields if,. The 
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by the spin Hamiltonian 

n 3 = -\ E JiASj - E = " 2 E j ^ s j - 2 E H & C 1 ) 

where Jjj = is presumed. Both quantum spin | operators and Ising spins Si(— 2s^) 
are used simultaneously in this work. Note that at this stage the bonds Jij are quite 
general. 

The assembly of spins is weakly coupled to a bath described by the Hamiltonian 
TL B via a spin bath interaction 

H SB = E B iSi = E \{Bfs- + Br 4) + B*s\ (2) 

i i 

where the operators Bi represent variables of the bath system. Thus the total 
Hamiltonian is given by 

h = h° + h sb with n° = n B + n s . (3) 

There is no need for an explicit specification of the bath Hamiltonian 7i B and the 
bath operators B^. As shown below it is just the absorptive part x'b^) °f the dynamic 
bath susceptibility 

POO 

Xb(uj) = x'bH + *x£M = -i ([Br,^ B *B+)) {B) e^dt , (4) 

Jo 

which enters the calculation and which is assumed to be known. In equation (^) (...)(£) 
represents the canonical thermal expectation value with respect to T~C B and L B denotes 
the Liouville operator defined by L B A = [H B , A] . For later use it is noted that x'b^) 
may be rewritten as 

= ± tt (e*** - 1) {B? 5(L B ± u) Bf) {B) (5) 

which can easily be shown or can be found in literature [Tj|. In writing equation (f|), it 
is assumed that the bath susceptibilities do not depend on the site %. This assumption 
is not essential and an extension to the general case is straightforward. 

Let tb and t$ be the relaxation times of the bath and the spin system respectively. 
Then it is natural to assume a fast relaxation to thermal equilibrium for the bath system 
compared to the spin system which implies 

r B < r s . (6) 
For the explicit determination of linear response quantities (see section 4.) the bath 
susceptibility Xb( u ) enters for ujts ~ 1. This implies ujtb 1 and 

Xb{ u ) ~ const u for u r^ 1 (7) 
can be used in the generic case for the explicit calculations of below. If the simple 
form Xb(u) = Xb(1 — ^I t b) is presumed, the constant of proportionality is given by 
const= Xb/tb, which may in principle be temperature dependent. This dependence is 
neglected in this work by assuming a slow variation on scale determined by the spin 
glass temperature. This assumption is to a certain extent arbitrary but can be justified 
for special cases. Such a case is the Korringa mechanism, where the bath and the Bi are 
identified as the conduction electrons and itinerant spin density operators respectively. 
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2.2. Mori formalism for the dynamic susceptibility 
The linear response of the magnetizations 

(Si)(t)-(Si) =Y,KMW expiut , t>0 (8) 

3 

due to small time dependent external fields h^ x Q(t) expiut is governed by the dynamic 
susceptibility matrix x{oj) which is given by the Kubo formula written in the Liouville 

space 



00 

Xi j (to) = P(S i \L{L + u + ir 1 }- 1 \S j ) , rj +0 (9) 

with A = A — (A). In the Liouville space the operators A of the state space are 
considered as vectors \A) with the temperature dependent Mori scalar product 

(A\B) = ^ / d\(A^e- xn Be xn ) = ^(A\l - e'^L^B) (10) 
P Jo P 

where (. . .) is the canonical thermal expectation value with respect to TC. Operators in 

the Liouville space like the Liouville operator 

L\A) = \[H,A]) (11) 

will be written in Roman letters H. 

Introducing the projection operator in Liouville space 

P = P^ t \S i )x^\S j \ with Xi j = P(Si\S j ) (12) 



and applying the standard Mori projection procedure [T4|, |15| leads to 

X-\u) = X- 1 -™p- 1 L-\u) (13) 

with 

Lij(u) = L'^iu) + iL'l^u) = i(Si\ LQ {QLQ + to + irj^QL] S 3 ) (14) 

and with Q = 1 — P. The matrices x an d L represent the static isothermal susceptibility 
matrix and the dynamic Onsager-Casimir matrix respectively. Note that the vectors | Si) 
which span the subspace P commute. Thus the frequency matrix vanishes in the present 
case. 

According to the standard approach [jTS for a sufficiently weak coupling Tt SB the 
Markovian approximation L(uS) w L'(0) can be applied in equation (0). Furthermore 
the leading order perturbation-theory expressions can be used for both x an d -^'(0) in 
(|13|). The static isothermal susceptibility matrix is of the order zero and approximated 
by 

Xij = (3 (Si Sj)( S ) = 13 ((Si Sj)(S) - mimj) (15) 

| The details of the Mori approach including the justifications and the general discussion of the 
approximations can be found in |p"5f . The notation of the present work widely agrees with |]l5| . Here 
both the Boltzmann constant and h are set equal to 1. 
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where rrii = (Si)rg\ and where (. . is the canonical thermal expectation value with 
respect to the Ising Hamiltonian ([[]). 

Setting L = L'(0) the lowest order of this matrix is the second order in the 
perturbation and is given by 

Lij = tt (Q°L SB Si | 8( Q°L°Q° ) Q°L SB ^) (0) = vr (L SB S l \ 8( L° ) L SB Sj\ 0) (16) 

where L SB and L° are the Liouville operators related to H SB and to 7i° respectively. 
The index (0) denotes the Mori product taken with Ti. alone and the projector Q° is 
defined with the latter Mori product. The projector Q° drops out, since Si commutes 
with H° and thus (Si\L B Sj) = i(Sj[H SB ,Sj]}^ = holds, where in addition the 
second equation of flTOl ) was used. Recalling Ju = 0, one finds that 

5( L° ) Bfsf = sf 5( L B ± 2H t ± 2X t ) B? (17) 

where the operator of the internal field Xi at the site i is given by 

Xi = JijSj • (18) 
j 

Using equation (|TT| ) we find = La Sij with 

7T / 

B O TT OV \ TD— i / 1 i C\ T~>- JT/T B 



(SiO)m = ( T 7 ^ ^7 ) = (Otanh/3(^ + ^)> ( s) (20) 



rt)^- 1 -^!)" 1 with — 5y ( 7loZr "I, L, (21) 



- ((1 - Bf8{h B - 2H t - 2X<) Bf + (1 + s i) BrS(L B + 2H i + 2X *) B t) -(19) 

iLi ' / (0) 

Let O be any operator in the unitary space of the spins not involving site z. Then 
TiiSiOe^-pHi -pXi 
Trt exp (-/?#* - Jx~) I (S) 
can be obtained ]r9[ . Applying this relation to equation ( |T^ ) and using (|5|) finally yields 
for the dynamic susceptibility matrix 

X%(2H i + 2X i . 
sinh /3(2Hi + 2XJ / (S) 
where the static isothermal susceptibility matrix x an d t ne internal field operators Xi 
are given by equation ( |l5l) and by equation ([18|) respectively. 

To complete the analytic investigations for the general Ising case, the local-field 
probability distribution functions 

P .(h) = ( S {h -Hi- Xi)) (s) = (S(h -H- JijSj))^) (22) 

are introduced which permits the La to be rewritten as 

™£m dh ■ (23) 

From the definition the relations 

J Pi{h) dh = l , J hPi(h) dh = H + Y^ Jij™j , J tanh(/?/i) P^h) dh = m, (24) 

immediately result, where the last relation is based on equation ( P0| ) with = 1. 

With the result ([H]) we have found in a very compact form of the dynamic 
susceptibility matrix x{ u ) ■ This result is not restricted to the SK model and holds for all 
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Ising models. According to the equations (pi) and ( p3|) the dynamic susceptibility x{ui) 
can be explicitly calculated provided that the static susceptibility x an d the internal 
magnetic field distribution functions Pi(h) are known. 

Knowledge of the linear dynamic susceptibility x( u ) implies knowledge of all the 
other response functions. Let us consider the linear relaxation functions <fr(t) which 
describe the linear response (Si)(t) = — Xij) hf for t > due to small changes 

of the external fields —h^ x Q(t). According to the general response theory &(t) is given 
by 

$ ij {t)=[3{S i \S j (t)) (25) 
and is approximated by 

**?Q = -PLx- 1 *(t) or by *(f) = X expi-x^fiLt) . (26) 

The remaining quantity of interest in linear response theory is the response function 
matrix which equals —<&(£). Due to this simple relation the response function matrix is 
not further considered in this work. 



The results ( p6|) can be compared in detail with the work of Szamel |12] which 
represents both the most recent and the closest treatment on the subject of this paper. 
Comparing equation ([56]) for the case H \ = with equation (13) of the first paper of 
reference [0 clearly shows that the values (1 — (Si tanh /3JQ}(5)) or according to relation 
( POD the values (1 — (tanh 2 fiX^rs)) are used for the La instead the correct values given 
in equation (pl|). Thus in general both approaches disagree. 

The work of Szamel and other former work [|11|] is based on the phenomenological 
Glauber dynamics. Microscopic derivation [T^, [17] of this master-equation (in the form 
of reference [ft^O leads to the transition rates 

-I'-^-wsi < 27 > 

for the Hi = case. In the work of Szamel the values for the rates (1 — Si tanh/3Xj)/2 
were used which can only be justified for the special, rather unrealistic case that the 
bath susceptibility Xb^) * s proportional to tanh(/3o;/2). Thus, for a general agreement 
with the present approach, the correct transition rates (27) have to be used in the 
master-equation treatment. 

Microscopically unjustified rates of the form of [I2|] are widely used to analyze the 
dynamics of Ising models for various physical questions. Provided the characteristic 
width of the distribution P(h) is of the order of the temperature T or larger than T a 
modified master-equation approach with the rates ( 27) will lead to significant changes. 
Note that for (standard) mean field treatments these effects are absent. In these cases 
P(h) is a 5-function and the modifications reduce to a factor which can be eliminated 
by scaling the time. For the spin glasses at low temperatures, however, the exact form 
of the transition rates is essential, as worked out in the following for the SK model. 
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3. Internal field distribution function for the SK spin glass 

For the rest of the paper we consider exclusively the special case of the SK model. In 
this case the bonds J^- are independent random variables with zero means and standard 
deviations iV~^. This scaling fixes the spin glass temperature to T = 1. The smallness 
and the randomness enter basically into the approach of this section where a tractable 
form for the field distribution Pi(x) of the SK model is deduced. 

This approach uses techniques similar to the derivations of the TAP equations 
H @, §. The spin Hamiltonian (D) is rewritten as H s = — (Hi + Xi)S{ + Hi where Hi 
describes aiV-1 spin system with the spin Si removed from the original spin system. 
Using the relation S(h — Hi — Xi) exp[/3(ifj + Xi)Si] = S(h — Hi — Xi) exp(j3hSi) and 
the Fourier representation of the ^-function and taking the trace over the site % 

= cosh^/i) r z t (k)e-^ h - H ^dk with Z i (k)=Tr l e^ J ^-^ (28) 

is obtained. The partition function of the original spin system is denoted by Z and Trj 
stands for the trace of the N — 1 spin system. 

The quantity Zi(k) represents a partition function of the N — 1 spin system in 
presence of additional imaginary fields ik(5~ l Jij. Due to the smallness of the J^- the 
cumulant expansion can be performed which yields to leading order in N^ 1 || [2|. || 

k 2 

In Zi(k) = In Zi(0) +ik^2 J ij™i ~ VqZ^ JijXjiJu (29) 

where rhj = f3~ 1 dlnZi(0)/dHj and Xji = drhj/dHi are the local magnetizations and 
the static susceptibilities of the N — 1 spin system respectively. Following again the 
former approaches || 0, [| the double sum in equation (|29p is replaced by the local 
static susceptibility 



1 x - 1 x drrij 



N ^— ' N ^ OH, 

j 3 J 

With the result (|29|) the integration in equation (|28|) can be performed yielding 



, (3 coshW r (3 X i\ / 0(h-H?")* y 
p i(h) = \ h ; ,„ rrP ffs ex P 1 ~ r ex P 1 o r ( 31 ) 



2n X i cosh(8H* ff ) P I 2 J I 2 X i 

where the factors independent of h are determined from the normalization condition 
1 = j Pi(h) dh and where the local effective field if? = Hi + Y2j Jij^j was introduced. 
The latter two relations of equation (|24]) yield 



H* ff = Hi + J2 J i3 m 5- m 5Xi and mi = tanii(j3H? ff ) (32) 
which is employed to rewrite equation (|3"T|) as function of h, rrii and X; 



//3(l-m?) r f /3(/i-/3 ^rtanhmi) 2 ^. 
P iW = U — — cosh(/?/i) exp | - — | exp | — j . (33) 
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Note that equation ([32]) are the well known TAP equations if the local static 
susceptibility \i is specified. In this work the modified version || of these equations is 
used which yields 

0(1 - ™?) 

Xi - 



Ly gilzgjj (34) 

ArZ^ i + r 2 /3 2 (l-m 2 ) 2 1 ; 



with 



r = for 1 - £ - ml) 2 > (35) 



N ^ 1 



As worked out in ||, for the stable states the modified TAP equations are equivalent 
to the original ones in the thermodynamic limit and differences result only for unstable 
states. For finite systems the complete temperature and field dependence of the solutions 
of the modified equations is known ||, in contrast to the original TAP equations. 

With the knowledge of these solutions the total internal field distribution function 
defined as 

i 

can be calculated according to equation (RTJ) or fl3"5|). In zero field above the spin glass 
temperature the distribution function reduces to 

P{h) = Pi(h) = cosh(/3/i) exp { - - + - } for T > 1 , H t = . (38) 

\/2tt I 2 J 



The latter result for the special case has already been obtained by Thomsen et al fl9 
However, to the best of the authors knowledge, the general results of this section have 
previously not been published. 



4. Discussion and numerical results 

4-1. Internal field distribution function 

The static properties of Ising models can entirely be formulated in terms of the internal 
field distribution function P(h) [fTS|, [T9H - Due to this general importance some discussion 
of P(h) for the SK model is presented. 

For zero magnetic field the temperature dependence of P(h) is shown in figure ([]]). 
Above the spin- glass temperature T = 1 the exact result ( |38"D is plotted. Below the 
spin glass temperature the figure is based on equation ( |3"3D and on the numerical results 
of 0. In the plot the local magnetizations m ; for the state of lowest free energy of a 
N=225 system (sample I) are used. On the scale of figure ([]]) both regimes T > 1 and 
T < 1 fit smoothly together. The distribution function P(h) bifurcates from a one- 
peak structure to a two-peak structure at the spin-glass temperature. With decreasing 
temperature the minimum located at h = becomes deeper and finally reaches the value 
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Figure 1. Temperature dependence of the internal field distribution function P(h) 
of the SK spin glass for temperatures T below and above the spin-glass temperature 
T = 1 for zero external field. 



P(h = 0) = for zero temperature within numeric precision. This behavior also applies 
to the metastable states. Moreover the variations for the different states are small and 
seem to become negligible in the thermodynamic limit. This is remarkable and is an 
indication for self averaging of P{h). These findings, as well as the general temperature 



dependence of P(h), are in agreement with Monte-Carlo simulations of references [19 . 

The discussion is completed by figure (§) where P(h) is shown for a system with a 
homogenous external field H = 0.5. No exact results are known for this case and thus the 
plot is based on numerical data everywhere. Again the data of sample I of reference [[J 
are used. Now the distribution P(h) is asymmetric but again a bifurcation to a two-peak 
structure is found when the spin glass regime is entered. This occurs at the Almeida- 
Thouless (AT) temperature [2C] which is determined from T 2 = N~ l 5^(1 — m 2 ) 2 and 
which leads to the numerical value of T = .577 for the sample under consideration. 
These results indicate that the spin glass regime is also characterized by a two-peak 
structure of P(h) for the case where a finite external field is present. 

Certainly in the numerical results finite-size effects are present. The most obvious 
feature is the asymmetry of P(h) in figure ([I]). Thus the investigations of the finite-size 
effects || are extended to the distribution function P(h) and averages over a few tens 
of independent samples are performed keeping N and T fixed. For the averages the 
asymmetry of P(h) for the case H = reduces both with increasing number of samples 
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Figure 2. Temperature dependence of internal field distribution function P(h) of 
the SK spin glass in presence of a homogenous external field H=0.5. The (numerical) 
AT temperature is T=.577 . 



and with increasing N . This represents some numerical evidence that the asymmetry 
is indeed an artifact due to the finite system size. 

In this context we recall that in zero magnetic field for each solution of the TAP 
equation (|32|) a further solution can be constructed trivially by changing the sign of all 
mj. According to equation (fTp) the distributions corresponding to these two solutions are 
given by P(h) and by P(—h) respectively. The means of these two distributions exhibit 
considerable smaller sample-to-sample variations than the individual P(h). This can in 
principle be used to construct improved approximations of P(h) for the thermodynamic 
limit. The TAP approach, however, does not use any averaging and thus we avoid this 
procedure. Moreover the asymmetry of figure (|l|) indicates the order of magnitude of 
the finite-size effects. 

4-2. Dynamic susceptibility 

We focus on the local dynamic susceptibility which is defined as 

XM = X 'i(u)+ix'i'(u) = N-'Trxiu) = N~ l Tr {x' 1 - iu^L)' 1 }- 1 (39) 

and which is a quantity of both theoretical and experimental interest. Employing the 
results of the last subsection and the approximation (|7|) (setting const=l which fixes 
the time scale) the Onsager coefficient La can be numerically determined according to 
equation (|23|). Furthermore using the well known expression for the static isothermal 
susceptibility 

X ^ = { [(3(1 - m*)]- 1 + xi] Sij - Jij (40) 
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with xi given by the equations ( |34| - p6D , all terms of equation (^) are explicitly known 
and the dynamic susceptibility matrix x{ u ) is obtained by numerical matrix inversion. 
From this the local susceptibility Xi{ u ) * s finally calculated. 

A further method exists to determine Xi{ u ) which is based on the theorem of Pastur 
PH [|. This theorem immediately leads to the identity (compare ||) 

X i(u) = N- 1 ^ - m])}- 1 + xi- tuJ^Lu)- 1 - ^(a;)}" 1 (41) 



and Xi( u ) is obtained as a solution of equation flflp satisfying ux"{u) > 0. In the 



paramagnetic regime and in absence of external magnetic fields this equation leads to 
the analytic result 



x,M = i±El^El£l + _ t r > 1>jri = (42) 

where the Onsager coefficient La = L{T) is temperature dependent and given by 

?2/0 , 1 f h exp(-/i 2 /2) 



L(T)^ex P (-^/2) y " («) 



using the equations ( P3|) and (p8|) . Note that for T ^> 1 equation ([42|) reduces to 
Xi = (3(1 — iojf3)~ l which implies a relaxation time proportional to /3, as is the case for 
the Korringa relaxation of single magnetic moments dissolved in a metal. Note further 
that the result (|42"D partially agrees with the work of Kinzel and Fischer |K| |3| who, 
however, have not obtained the temperature dependence of L(T). Near the spin glass 
temperature both approaches agree and yield Xi{w) = 1 + at T = 1. 

With reference to the definition of Y p|. || 

Xi(oj) | w _±o= xi ± (44) 
results from equation (fi~l~|). Thus the quantities xi and T, which enter formally in 



derivation the modified TAP equations, are related to the real and the imaginary part 
of the isolated susceptibility xi( u ~> 0)- Recalling that there are frequently situations 
in physics where the isolated and the isothermal susceptibility differ (for some general 
discussion see |15|) this result gives some additional insight into approach leading to the 
modified TAP equations [|, [J. 

In figure (0) both the analytical and the numerical results for xi are presented for 
the zero external field case. For the numerical parts again the lowest free energy state 
of sample I from reference || with iV = 225 is used. The two different methods lead to 
similar results but deviate increasingly for small frequencies u. It is remarkable that the 
results of figure (§) show some overall similarity to real, experimental data [§, |J. One of 
these features is the frequency dependence of the cusp temperature of x'i(T, u = const). 
In the present approach the shift results simply from the temperature dependence of 
the Onsager coefficients. 

§ Actually a slight generalization of this theorem is needed to incorporate the term iuifl -1 L^ 1 . Such 
a generalization can be easily be added to the method of Bray and Moore |22l . 
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Figure 3. Temperature dependence of the real part, and the imaginary part, 

Xi'(to), of the local dynamic susceptibility xi( w ) m zero external field for different 
frequencies. The data points are calculated by numerical inversion of equation (fiol ) for 
a system of N = 225 spins (compare text). The lines represent the results calculated 
from equation ( fill) y i a the theorem of Pastur. For T < 1 they are again based on the 
numerical data of the N — 225 spin system. For T > 1 the lines give the analytic result 
of equation ( [l2| ) in the thermodynamic limit. The dashed line, the dotted line and the 
full lines correspond to u> — 0.001, to co = 0.01 and to w = 0.1,1.0,5.0 respectively. 
For ui = 0.001 some data points exist which are outside the plot-range of Xz( w )- 

The deviations appearing in figure (|3|) for small frequencies are caused by the 
finite-size of the system. According to equation fl2~T|) L is positive definite. Thus 
the matrix A = {(3L)^ is well defined and x _1 (u;) can be rewritten as X _1 (^) = 
A _1 (Ax _1 A — iuj^Ar 1 . Use of the diagonal representation of Ax -1 A leads to 

Xi(w) = N- 1 ^2(u a \A 2 \u a ) (6 a - iuj)' 1 where Ax^A\u a ) = 5 a \u a ) (45) 

a 

which shows that Xl( u ) can De written as a superposition of Lorenzian functions with 
relaxation rates S a given by the eigenvalues of Ax" 1 A . In the spin glass regime the 
minimum value 5 m i n of all rates 5 a is small (the numerical value for the data used in 
figure (§) is of the order 0.03) for finite N and tends to zero in the thermodynamic limit 
H |, p2| . For frequencies uj < 5 m i n the susceptibility X;(u;) sensitively depends on the 
small eigenvalues according to equation (|45|) and thus finite-size effects mainly show up 
for small uj for the data obtained by direct matrix inversion. According to figure (|3|) 
the deviations are moderate for those results which use the theorem of Pastur. This 
indicates that the use of this theorem for finite iV smoothes partly out the finite-size 
effects. 
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Figure 4. Time t dependence of the reduced local relaxation function E(i) = 
(0) in a double log scale for different temperatures T in zero external field. 
The data points correspond to the numerical results of a sample with N — 225 spins. 
The lines represent the results in the thermodynamic limit given by equation ( [f7| ) . The 
full, the dashed and the dotted line correspond to T = 1.0, 1.2 , and 2.0 respectively. 



4-3. Relaxation function 

The local relaxation function can be written as 

$,(*) = N^Tr&it) = N- iy 52(u a \A 2 \u a ) 5 a - 1 exp(-5 a t) (46) 



where equation (2B) and the eigenvalue equation Q4"5|) and were used. Again the 
superposition of the contributions resulting from the different rates 5 a (or relaxation 
times 5~ l ) can be identified. The local relaxation function satisfies 3>j(t) < $z(0) = Xz- 
and thus is well behaved in the thermodynamic limit even for the case when 5 mi „ tends 
to zero. For vanishing external fields $/(£) can be given analytically 

Ht) = lf d X f3 ~ X * 2x exp [ - (1 + (3 2 - 20x)L(T)\ T> 1,^ = (47) 

in the paramagnetic regime [l|]. As already noted, the temperature dependence of L(T) 
differs from Jj], |1^] and is given by equation (|43|). 

In figure (^) the numerical results for the local relaxation function $j(t) are plotted 
(again for the state with lowest free energy of sample I of reference 0) together the 
analytic results fl4"7p. The plot clearly exhibits the slowing down at the spin glass 
temperature and the presence of the slow dynamics in the spin glass regime. Comparing 
the numerical and the analytical results above the spin glass temperature the deviations 
increase for the long time behavior when approaching the spin glass regime. As discussed 
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already above, these finite-size effects result from the finite value of 5 m i n and limits the 
numerical results to the region t <C 5 m i n ~ 30 near and below the spin glass temperature. 

5. Conclusions 

Two questions are studied from a more general point of view in the present work. First 
of all the linear response theory for an arbitrary Ising model which is microscopically 
and weakly coupled to a bath is investigated using the theory of Mori and applying 
the standard approximations. This approach seems to be natural and conservative and 
could have been carried out some decades ago. Nevertheless the results are remarkable, 
as the entire dynamical response of any Ising model is completely determined by the 
internal field distribution function, by the static isothermal susceptibility and - as the 
only characteristic feature of bath - by a bath dynamic susceptibility. Due to the simple 
structure this part of the present work may be of interest for other models than the SK 
model. 

The relation to former approaches has been worked out in some detail at the end 
of section 2. It is the use of microscopically unjustified transition rates in these master- 
equations which causes the differences to the present investigation. 

The second question is exclusively related to the SK model and represents the 
study of the internal field distribution function showing how this function is related to 
the solutions of the TAP equations. The relationship obtained opens the possibility 
of explicitly calculating further quantities of interest for the SK model, such as the 
inelastic-neutron-scattering cross section 



A further tool to obtain the final results for the dynamic linear response functions of 
the SK model is the explicit knowledge of the solutions of the modified TAP equations. 
This again demonstrates the importance of the modified TAP approach. 

The present work is limited to the dynamics in linear approximation near the 
thermodynamic equilibrium. It is, however, well known that in the physics of spin glasses 
nonlinear dynamical effects and the linear response in out of equilibrium situations are of 
great importance ||. Thus an extension of the present microscopic approach to nonlinear 
dynamics would be of great use to treat these effects on a well-founded bases. Such an 
approach based on the theory of Robertson or on the theory of Nakajima-Zwanzig (see 
e.g |15|| ) is in progress and will be published separately fl6| . 
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